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In this paper we explore a possibility that all transport turbulent models are contained 
in a coarse-grained kinetic equation. Building on a recent work by H.Chen et al (2004), 
we account for fluctuations of a single -point probability density in turbulence, by in- 
troducing a "two- level" ( c, v)-phase-space, separating microscopic (c' = c m i cro = c — v) 
and hydrodynamic (v' = v — V) modes. Unlike traditional kinetic theories, with hydro- 
dynamic approximations derived in terms of small deviations from thermodynamic equi- 
librium, the theory developed in this work, is based on a far- from -equilibrium isotropic 
and homogeneous turbulence as an unperturbed state. The expansion in dimensionlcss 
rate of strain leads to a new class of turbulent models, including the well-known JC — £, 
Reynolds stress and all possible nonlinear models. The role of interaction of the fluxes 
in physical space with the energy flux across the scales, not present in standard model- 
ing, is demonstrated on example of turbulent channel flow. To close the system, neither 
equation for turbulent kinetic energy nor information on pressure-velocity correlations, 
contained in the derived coarse-grained kinetic equation, are needed. 



1. Introduction 

In this paper we revisit an old and all-important problem of turbulent modeling. The 
problem, first formulated in terms of images borrowed from kinetic theory by Prandtl 
(1925), was later developed by Kolmogorov (1942), Launder and Spaulding (1974) and 
many of their followers. The success of these semi-empirical approaches can hardly be 
overestimated: with development of powerful computers turbulent modeling became an 
important part of scientific and engineering design process. Still, these models, involving 
low -order time - derivatives and low powers of dimensionless rate of strain, are less effec- 
tive in describing strongly sheared and rapidly distorted flows where the dimensionless 
rate of strain is not too small. To improve performance, various non-linear models, pio- 
neered by Speziale (1987), accounting for the next order in powers of the rate of strain, 
have been proposed. Since the second-order trancation of a lacking-small -parameter ex- 
pansion is problematic, these models had a mixed success. (The role of the large-eddy 
simulations (LES) in description of complex flows will be discussed in Conclusions. ) 
The attempts to systematically derive turbulence transport models from the Navier- 
Stokes equations (NS), based on renormalized Wyld's perturbation expansions (RNG 
(Yakhot/ Smith et al (1992), Rubinstein and Barton (1990)), double expansion (Yakhot et 
al (1992)), DIA (Yoshizawa (1987)), were relatively successful and the resulting equations 
found their place in engineering. It became clear soon that resummation of the series was 
impractical and the approach was restricted to the low- orders. 

It has been shown recently by H. Chen et al (2004) that the turbulence models, derived by 
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application of the Chapman-Enskog expansion (CE) to a model kinetic (Boltzmann-BGK) 
equation written for hydrodynamic modes, are similar to those previously obtained from 
the Wyld (1961) expansion, applied directly to the Navier-Stokes equations. The main 
result of Chen et al (2004) can be stated as follows: the perturbation theory, leading 
from kinetic equation to the Navier-Stokes (NS) equations, generates the well-known 
turbulence models if, instead of the relaxation time t — v a /8 — const of kinetic theory, 
one writes t t K/£ and represents the temperature 9 as 9 = K. Here K and £ are 
the mean turbulent kinetic energy and dissipation rates, respectively. Thus, the formal 
superficial similarity between the two perturbation expansions has been established. This 
result means that, in principle, the models of an arbitrary non-linearity and complexity, 
which we even cannot explicitly write down, are contained in a simple kinetic equation. 
In this work, we would like to reformulate the procedure developed in Chen (2004) 
by taking into account both microscopic and hydrodynamic modes. The coarse-grained 
kinetic equation for hydrodynamic modes is then obtained by integrating out the fast 
microscopic variables. To achieve this goal, we had to substantially redefine the expansion 
procedure and, as a result, obtain a qualitatively new class of turbulence models. 

1.1. The Boltzmann -BGK equation. 
Kinetic theory for a low-density gas is based on the Boltzmann equation 

|[ + cV/ = C(/) (1.1) 
for a single particle distribution function /(c,x, t), where collision integral 



c(f) = / w'(f'f[ - f h) drw'dr; = / \c rel \(ff[ - f f 1 )dad 3 Pl = i in - w (i.2) 



If the system Hamiltonin is known, this equation can be derived directly from the Liuville 
theorem. Here the element of the phase -space volume dT = dxdp and \c re i\ = |ci — c| is 
the relative velocity of colliding particles. If the intermolecular distance is A >> a where 
a is the length-scale of intermolecular interaction, the strait particle trajectories between 
collisions are assumed in deriving expression (1.2). The probability density functions in 
(1.2) are /' = f'(r,c',t) and /{ = f'(r, c'-y,t). Since in accord with the theory of elastic 
collisions, the velocities c', c' x are expressed in terms of c, Ci, the integration in (1.2) is 
carried out over the phase space Li only. The relaxation -time approximation (RTA): 

C = - J - J — (1.3) 

T 

where f eq is the distribution function in thermodynamic equilibrium and r is a properly 
chosen relaxation time, is an often used anztaz, mentioned in Landau and Lifshitz (1981) 
as a "rough estimate of the ( Boltzmann) collision integral" (1.2). Still, recent implemen- 
tations of this approximation in the so called "Lattice Boltzmann" numerical codes, led 
to a remarkable success in simulating a wide variety of extremely complex fluid flows. As 
of today, the approach has been tested on basically all examples of classic laminar flows 
(Benzi (1992), Chen (1998), Succi (2001)). The robustness, speed and simplicity of the 
method, made it an attractive tool for both theoretical investigations and engineering 
design. 

Even more spectacular, and somewhat unexpected, is success enjoyed by the Lattice 
Boltzmann method in simulating complex turbulent flows (H. Chen (2003)). In this 
application, the relaxation time r is expressed in terms of hydrodynamic observables: 
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t w /C/£ , where turbulent kinetic energy is /C = v^. ms /2 and the mean dissipation rate 
is £ = z^(|^) 2 . In this case, in the first order of the Chapman-Enskog (CE) expansion 
one obtains the well-known JC — £ model with the non-linear JC — £ model appearing 
in the next order (Chen et al (2004)). In addition to speed and simplicity, the most 
attractive feature of this type of turbulence modeling is that it does not require theory 
of pressure-velocity correlations which, in complex flows, is an extremely difficult, still 
unsolved, problem. 

Typically, derivation of hydrodynamic Navier-Stokes equations from kinetic theory is 
based on a few assumptions. 1. If spatial gradients of mean velocity, density and temper- 
ature are equal to zero, the gas is assumed in thermodynamic equilibrium; 2. Expansion 
in powers of dimensionless rate of strain r] = rjij — rSij where Sij = \ + g^r) leads 
to transport equations describing various physical phenomena. 

The equations, derived this way can be used to fully describe large Reynolds number 
turbulent flows requiring at least N w Re 3 number of degrees of freedom. In the flows 
of practical importance this number is huge and therefore, the value of a coarse graining 
procedure leading to turbulence models can hardly be overestimated . 
It has been shown (Yakhot et al (1992) ) that the systematic derivation of turbulence 
models directly from the Navier-Stokes equations, can be formulated in terms of Wyld's 
(1961) diagramamtic expansion in two dimensionless parameters: the Reynolds number 
Re = u rms L/vQ and, familiar from kinetic theory, dimensionless rate of strain 77 = ttS^. 
The first expansion is responsible for description of isotropic and homogeneous turbulence 
with the non-zero energy flux across the scales and the second - for the non-zero spatial 
fluxes and flow structures. Due to proliferation of tensorial indices in the high-order 
contributions, the resummation of the expansion is an extremely difficult task and it has 
been conjectured by Polyakov (2001) in the middle of seventies that it may result in a 
kinetic equation containing all terms of Wyld's scries. If this is true, a turbulence model 
of an arbitrary non-linearity and complexity may be contained in a relatively simple 
kinetic equation. The most important and interesting feature of this system is that, unlike 
equilibrium, the r\ = state of isotropic and homogeneous turbulence is not flux- free but 
involves a large 0(1) energy flux across the scales. Thus, to derive turbulence models, it is 
desirable to develop a kinetic approach not based on the equilibrium, flux-free, zero-order 
state. 

In this paper we, building upon a remarkable work of H. Chen et al (2004), consider both 
turbulent (hydrodynamic) and thermal (microscopic) velocity fluctuations as governed by 
the Boltzmann -BGK kinetic equation. Since microscopic and hydrodynamic fluctuations 
occupy their respective fractions of the phase-space, this equation is formally defined on 
a somewhat enlarged phase-space. The coarse-graining, eliminates the small-scale fast 
microscopic fluctuations and restores the six-dimensional space, leading to the kinetic 
equation for the large-scale (turbulent) velocity component only. In this case, the zero- 
order, zero-mean-spatial-gradient state is not an equilibrium but that of isotropic and 
homogeneous turbulence characterized by a finite energy flux across the scales. The 
Chapman-Enskog expansion applied to this state leads to a novel set of turbulence models 
of arbitrary nonlinearity. 

This paper is organized as follows. In the next Section we, for the sake of clarity and 
continuity, describe traditional derivation of hydrodynamic approximations from kinetic 
equation and briefly outline the way turbulence models appear from the NS equations. 
In Section 3, which is most important for the present development, we introduce the 
non-equilibrium kinetic equation and the zero-order pdf for isotropic and homogeneous 
turbulence, both defined on an enlarged phase-space. In Section 4, the finite energy flux 
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is defined as a dynamic constraint and it is shown that the introduced in Section 3 non- 
equilibrium pdf is indeed a solution to non-equilibrium equation which conserves the total 
kinetic energy of a system. Finally, in Section 5, the coarse-grained kinetic equation for 
hydrodynamic modes is derived. It is shown in Sections 6 and 7, that the CE expansion, 
leads to transport equations for both turbulent velocity field and kinetic energy In Section 
8, choosing a proper relaxation time, we show that our kinetic equation contains the well 
known K-E , Reynolds stress and non-linear models. In addition, this simple equation 
contains turbulent models of an arbitrary non-linearity and complexity. Some new effects, 
originating from interaction of spatial fluxes with the energy flux (cascade) arc identified. 
Summary and conclusions are presented in Section 9. 

1.2. Equilibrium. 

In equilibrium, with all temporal and spatial gradients equal to zero, the left side of 
kinetic equation (1.1), (1.2) is equal to zero and the remaining equation C = has a 
solution: 

E(r) 

f eq = constant x e — (1.4) 

where the energy E(T) is conserved in the process of collision, i.e. E' + E[ = E + E\ and 
the temperature 9 = const. For a gas moving with the constant velocity V, the solution 
is: 

E(r)- P -v 

f e i = const x e 5 (1.5) 
This solution is clear if we introduce internal energy e in and write 

m < 

E(T) = e in + — (1.6) 
In a frame of reference moving with velocity V, we have 

f eq — constant x e 5 e 55 (1.7) 

In deriving transport approximation, one expands the probability density powers the 
large-scale field (LL) V w V + §^~ x i- The constant component Vq = const can be 
removed by the Galileo transformation and only derivatives of the large scale field V are 
relevant. 



2. Hydrodynamic approximation 

In this Section we, for the sake of clarity and continuity, briefly describe the main steps 
leading from the Boltzmann equation to hydrodynamic approximations. 
Classic derivation of transport approximations starts with the Boltzmann equation (1.1), 
(1.2) or (1.3) and equilibrium pdf 

r = ^exp(- (c - V o ( ;^ ))2 ) (2.1) 

In this Section we denote velocity of the frame of reference v(x, t) which is a physical 
(not phase-space variable). In what follows we restore the notation (V) used above. 
This should not lead to any difficulty. The equilibrium pdf is obtained from equation 
C(f eq ) — 0, which is valid locally, so that the hydrodynamic field v(x, t) can be a slow 
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varying function of time and space coordinates. In addition, the conservation laws and 
the symmetries of collision integral give: 

J Cdc = J C(f)cdc = (2.2) 

Introducing the fluid density p = J f(v)dv and integrating (1.1)-(1.2) over c leads to the 
continuity equation: 

^ + V-pv = (2.3) 
Multiplying (1.1)-(1.2) by c with subsequent integration over the phase space yields : 

+ dipvjVi + di<Jij = (2.4) 



dpvj 



dt 



where the stress tensor is : ov, = p(a — Vi)(cj — Vj). In thcrmodymanic equilibrium, the 
stress-tensor is evaluated readily with the distribution function f eq from (2.1) giving 
°Xj |e<7 = pOSij, so that combining this with the ideal gas equation, gives: 

d 1 

-g-+VidiVj--Vp = (2.5) 
at p 

In what follows A\ eq = J Af eq dc. Defining 9 = i(cj — Vi) 2 , we obtain: 
89 



— + Vipvi9 + -Wip(ci - Vi)(cj - vj) 2 + ^p( c i - Vi)(cj - Vj)Sij = (2.6) 

On the equilibrium pdf f eq , = 95ij and the third-order in fluctuations contribution 
to (2.6) is equal to zero. 

2.1. Turbulence modeling. 

The nature of the large-scale hydrodynamic field v in (2.4) has not been specified. The 
problem of turbulence modeling is formulated as follows: let the velocity field v = V + v' 
where v' = and v = V. Our goal is to derive the equation for the large-scale slow field 
V. The above decomposition gives: 

dV 



+ V- VV + v'-Vv' + -V- ( r = (2.7) 
at p 



or 



dV 1 

— h V • VV H — V • (a + a v ) — v'V • v' = (2i 

at p 



where 



<Tij,v = (y-V) i (y-V) j (2.9) 

Now we need a theory or at least set of rules to express both stress a and the Reynolds 
stress v' • Vv' in terms of the large-scale field V and its derivatives. This program has 
been systematically developed in Yakhot et al (1992), Smith et al (1992). The method 
consists of a few steps: 1. In the first order, the stress tensor aij w 2voSij. This leads to the 
Navier-Stokes equations (NS). 2. Add a statistically isotropic and homogeneous stirring 
force to the right side of the NS. This force generates an 0(1) energy flux V = f • v' = £ 
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Figure 1. Scale and velocity field intervals in turbulence: 1. Energy range : r > L (V ); 
2. Inertial range: t\k < r < L (v'=v-V). PDFs: r « L (Gausss); r <C L (Multifractal); r < X 
(Microscopic) (c'=c-v ). 

across scales. In the case of a free turbulence with the total kinetic energy K, prepared 
as an initial condition, the energy balance reads d t K — —£. 3. Use Wyld's (1961) 
expansion to eliminate small-scale velocity fluctuations. It has been shown (Yakhot et al 
(1992)) that development of turbulent models involves a double expansion in powers of 
two dimensionless parameters: the Reynolds number Re = v' rms L/v and rate of strain 
t]ij = rj = rSij where tt is the relaxation time of turbulent fluctuation. The expansion in 
Re leads to the fixed point i?e* w v' rms L/vT = const, corresponding to the zero- spatial 
-gradient state of isotropic and homogeneous turbulence which serves a a ground state 
for the additional expansion in powers r\. The first expansion accounts for the non-zero 
energy flux of isotropic and homogeneous turbulence while the second - gives rise to 
the not zero mean velocity field and various fluxes. In the first order of the ^-expansion 
one derives a variant of JC — £ model, in the second -the Reynolds stress model and 
non-linear models etc. It is important that the second order 0(iy^/KS ■ S) terms of the 
non-linear turbulent models, are well known from kinetic theory as next-order in the 
Chapman-Enskog expansion corrections to the Navier-Stokes equations (Landau-Lifshitz 
(1981) ). This fact establishes similarities between Wyld's expansion of the NS equations 
and the CE expansion of kinetic theory. While, one can formally write each term of the 
^-expansion, resummation of an infinite series is totally out of reach. 



3. Formulation in terms of kinetic equation. 

In turbulent flows, the velocity field is characterized by three correlation lengths (Lan- 
dau and Lifshitz (1959), Frisch (1996)): 1. microscopic mean -free path r w A w 
c s t, where c s is the speed of sound and r is the microscopic relaxation time; 2. hy- 
drodynamic inertial interval t/k ~ LRe~^ < r < L, where both large-scale forcing 
and viscous effects can be neglected and 3. nonuniversal, geometry-dependent, energy- 
production range r > L. We denote three velocity fields corresponding to these intervals 
as c' = c m i cro = c v; v' = v — V; V, respectively. (See Fig.l). Locally, at the scales 
r w A, the system is in thermodynamic equilibrium so that the c m i cro -field obeys Maxwell 
(Gibbs) statistics. 

Hydrodynamic approximations follow elimination (small-scale averaging) of modes from 
the interval r < A w c s r, where r is the relaxation time characterizing dynamics of 
small perturbations from thermodynamic equilibrium. After that, the turbulence models 
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are derived by averaging the Navier-Stokes equations in the interval t]k < r < L. We 
do not doubt the ability of the Navier-Stokes equations fully describe turbulent flow 
of arbitrarily large Reynolds number. In this case, since the number of excited modes 
is at least 0(i?e 9 / 4 ), the computational cost of these direct simulations (DNS) is very 
high. Our concern here is with derivation of the coarse-grained equations (turbulence 
models) for partial representation of turbulence in the interval r s=s L, which is hard to 
obtain from hydrodynamics. In other words, in this work we are interested in the mode- 
elimination from the interval < r < L applied directly to the Boltzmann equation. 
If successful, this way we will obtain the coarse-grained kinetic equation, containing, as 
a hydrodynamic approximation, all possible turbulence models, which are very hard to 
derive from the Navier-Stokes equation. In particular, we are interested in equations 
of motion for the coarse-grained or filtered velocity field which includes only a small 
fraction of hydrodynamic modes corresponding to the top of inertial range, i.e. r w L. 
Whatever the filtering procedure, the coarse-grained fields are obtained from the full 
fields by averaging over flow volumes of linear dimension r « A < L, where A is the 
filtering scale or over the modes v(k) with k > 1/A in a Fourier space. Also, one can 
define a coarse - grained field as: 



where the integration is carried out over the volume 1Z w A 3 centered about the point x 
or 



where F is a filter function. We see that the coarse-grained velocity field is a multipoint 
construction and evaluation of high-order moments is not a simple task. Moreover, spatial 
velocity derivatives are defined on a cut-off as: 



and, if A <C L, their moments, called structure functions, obey multifractal statistics. The 
inertial range multifractality (anomalous scaling ) means that the moments of different 
orders are independent of each other. As a result, the full statistical description of velocity 
field involves a large, infinite in the limit Re — > oo, number of independent parameters. 
This makes an accurate coarse-graining with the filter (mesh) size A << L basically 
impossible. 

The situation changes dramatically in the limit A — > L. It is a well-established, though 
not well-understood, fact that at the scales r — A s=a L, both the turbulent velocity fields 
v' = v — V and Sv' = v'(x + A) — v'(x) obey close- to- gaussian statistics. Moreover, 
it was demonstrated by Schumacher et al (2007) that in the low-Reynolds number tur- 
bulence, the velocity derivatives, too, are Gaussian random variables. If this is so, the 
entire large-scale field can be accurately described in terms of the second-order moment 
only. Our goal is to derive a coarse - grained equation for the modes v' fluctuating on 
the scales r w L. In the standard turbulence modeling, this is operationally achieved by 
introducing a large "effective viscosity" vt ~ u rms L, where 2K = uf ms w ^CkL 2 , into 
the Navier-Stokes equations. Due to a particular choice of viscosity vt S> vq, the small- 
scale modes (r <C L) are overdamped and the resulting hydrodynamic equations describe 
large-scale properties (r L) of the flow only. The defined this way equation is called " 
transport model" as as opposed to the LES models based on an inertial range cut -off 




v(x) 




d x v = (v(x + A x ) - v(x))/A : 



x 
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r = A«L. Since the dominant contribution to turbulent kinetic energy is contained in 
the large -scale fluctuations, we write the zero-order nonequillibrium probability density 
function in the laboratory reference frame: 

In this expression, both c and v- fields are phase-space variables independent of time 
and space. Thus, we have temporarily enlarged the phase space from six to nine di- 
mensions. The flow described by the probability density function (3.1) can be perceived 
as an ensemble of laminar flows each in thermodynamic equilibrium with the fluctuat- 
ing equilibrium Maxwell-Gibbs probability density. It is clear that: / f°dvdc = p and 
c| = / cf°dcdv = V and v| = V. In the laminar limit JC — > 

and integrating over v, we recover the expression (2.1). The pdf (3.1), resembling that 
of a mixture of gases with different temperatures and JC, describes an essentially non- 
equilibrium state. It will be shown below that (3.1) is a solution to a "nonequilibrium" 
kinetic equation and this temperature difference is related to the energy flux from hydro- 
dynamic modes to the microscopic ones. 
From the above definitions: 



V = v| = J v/°(c,v)dcdv 
£°(v)/°(v)= J c/°(c,v) ( ic = v/ (v) 



£ = i(v- V) 2 | = i/(v- V) 2 /°(v. c 1,/cr/v 



(3-2) 



In what follows, we assume the probability density 



/(c,v,x,t) = p/(^,^=X,t) (3.3) 



so that locally (c — v) 2 = d9 and (V — v) 2 = dfC 

The pdf of hydrodynamic velocity field /(v) is obtained from (3.3) by integrating out 
the microscopic variable c, i.e. 

/(v) = /"/(<:, v,x,t)dc (3.4) 

Our goal is to derive a kinetic equation for /(v) describing turbulent fluctuations only. 
To achieve this goal, we first increase the original six- dimensional phase space of the 
previous section (c,x) to the nine- dimensional one (c,v,x) and, eliminating thermo- 
dynamic fluctuations, project the nine- dimensional problem to the six- dimensional one 
(v, x). The advantage of this procedure from the turbulence modeling viewpoint is that 
while original problem involved N\ rts fl/X 3 number of particles, the coarse-grained one 
- only Nl w Q/L 3 , where the integral scale L >• A is the large scale correlation length 
of turbulence. 

Since the fast c m j cro = c — v -field varies on the time scale r w A/c s , where c s stands 
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for the speed of sound, (c — v) • v — 0, we can consider these components statistically 
ortogonal and introduce the continuity equation in the enlarged 9d-phase- space. In this 
case the Botzmann equation (Landau/Lifshitz (1981)): 



where /(/) is not yet specified 'non-equilibrium collision integral' with /°(v, c, x, t) ob- 
tained from equation: 



with V — const (Sij = 0) . Thus, in the theory developed below, we abandon the equilibrium 
pdf (2.1) for the the zero- spatial- gradient state of a gas in favor of the nonequilibrium 
expression (3.1). The detailed discussion of this point will be given below. Here we 
just mention that due to the non-zero energy flux across the scales, the kinetic energy 
of the zero-spatial-gradient (77 = 0) turbulent flow decays leading to increase of the gas 
temperature 9. 

4. Finite energy flux as a dynamic constraint. 

If a turbulent flow is supported by a large-scale energy input per second V, by nonlin- 
earity, this energy is transfered to the smallest scales where it is dissipated. It is usually 
assumed that the nature of the dissipation mechanism is unimportant. On the inertial- 
range scales, where both energy pumping and dissipation are negligibly small, the dynam- 
ics are characterized by the mean energy flux J = V = £ , where £ is the mean dissipation 
rate. In case of isotropic and homogeneous turbulence, prepared at initial instant of time 
with kinetic energy K(0), the energy K(t) decays to zero. Since in this case, all spatial 
gradients are equal to zero, no mean spatial fluxes are involved in the relaxation pro- 
cess. Unlike equilibrium gas, this strongly non- equilibrium system is characterized by a 
non-zero constant energy flux across the scales . 

It has been shown in a series of extraordinary works by Zakharov et al (1975), (1992), 
(1984) that, in case of turbulence, the Euler equations, with velocity field written in 
terms of Clebsch variables, can be represented as a kinetic equation for waves (particles) 
with the collision integral Iq describing the non-linear interactions. It has been shown 
that, in addition to equilibrium solutions, the equation Ic{n(k)) = for the occupation 
numbers n(k), playing the role of the pdf in the /c-space, has solutions corresponding to 
the non-zero constant fluxes in the wave-number space. The direction of these fluxes is 
exclusively determined by physics of the problem and peculiarities of non-linear interac- 
tions. For example, the mean field approximation developed in Yakhot (1992), led to the 
Kolomogorov energy spectrum as one of Zakharov's constant-flux solutions. Thus, unlike 
thermodynamic equilibrium, the kinetic equation for a zero -mean-gradient turbulent 
flow must satisfy the flux (energy balance) constraints: 



^-+c-V/ = J(/(c,v)) 



(3.5) 



-^=J(/°(c,v))^0 



(3-6) 



V 



£ 



(4.1) 



valid in the statistically steady state or 



dIC 
dt 




(4.2) 



- in case of decaying turbulence. 
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The relations (4.2) and (4.3), reflecting the energy balance in a flow, are independent 
upon the nature of the dissipation mechanism. We split collision integral /(/) into two 
components: Ihi(f) responsible for the non-zero energy flux of isotropic and homogeneous 
turbulence and I 1 (f) - for the spatial fluxes of a general turbulent flow. Since, as was 
pointed out above, the mean energy flux is a large-scale property, no differential operator 
can enter the expression. Thus we have: /(/) = I hi + 1 1 and 

/ r (v-V) 2 -d/C AC (v-c) 2 -d0.£ ,. , , in . 

4.1. Isotropic and homogeneous turbulence. 

Now we will show that the zero-order pdf (3.1) is a solution to kinetic equation describing 
isotropic and homogeneous turbulence. Since in this case, all spatial derivatives are equal 
to zero, the kinetic equation is: 

df° 

■w = I » (44) 

and c = v = 0, c 2 = dO and v 2 = dJC. Substituting (3.1) into (3.5)-(3.6) gives: 
!/( v - c ) 2 , d0 2., 1 ,v 2 ,..d/C 2„ 

/-T^ - rf % - = Jz { lc - d){ m + d £) (45) 

This equation, valid for all values of independent parameters c and v, can be correct only 
if: 

If- 

Thus, we conclude that the zero-order pdf (3.1) is indeed a solution to the kinetic equation 
(3.5), (3.6) for isotropic and homogeneous turbulence, subject to the energy constraints 
(4-1), (4-2). One can also see that according to (4-6), (4-7): |(/C + 0) = K + T = const, 
which reflects conservation of the total kinetic energy in decaying isotropic and homoge- 
neous turbulence. 



5. Coarse-graining. 

Introducing conditional probability density Q(c|v) we can write: 



so that 



/(v,c,x,t) = /(v,x,t)Q(c|v) (5.1) 

/(v, x, t) = j /(v, x, t)Q(c|v)dc (5.2) 
c a Q(c\v)dc = f a (v) + v a (5.3) 



and 



J dvCa(y)f( 



v) = V a (x) (5.4) 
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Below we develop an approach resembling the one proposed for the probability densities 
of the passive scalar in a random velocity field in Sinai, Yakhot (1989). Integrating the 
nine-dimensional Boltzmann equation (3.5), (4.3) over the relatively fast- field c gives: 

^ + V • v/(v) + V • £(v)/(v) = 1(f) + i 1 (5.5) 

where 



and the collision integral 



,i_ /~/°(v) 

TT 

is formally written with the yet unspecified relaxation time tt- The coarse -grained zero- 
order pdf is: 

/°(v) = T^Tire-^ (5-7) 

(27T/C) 2 

5.1. Model for conditional expectation value. 

The expression for conditional expectation value £»(v) is a subject to four dynamic con- 
straints. 

1. It must be invariant under Galileo transformation and not violate the Galileo invari- 
ance of equation (5.5). This is the reason for the v a term in the right side of (5.3). 

2. Integrating (5.5) over v must lead to the continuity equation. This gives: 

V • J £(v)/(v)dv = (5.8) 

3. Multiplying (5.5) by v and comparing the result with (2.6) yields: 

V 4 J Vjti(v) f°(v)dv = \7 lP (c t -v l )(c J - Vj )\ = V jP 6 (5.9) 

4. In the zeroth order of expansion in powers of dimensionless spatial gradients (Chapman- 
Enskog expansion), the equation (5.5) must generate (2.5) with the pressure term con- 
taining information about the integrated out microscopic variables. In a simplest case of 
ideal gas, p' = p9. 

A relation satisfying the above constraints is: 



£(v) = • W(0) (5.10) 



where 



s ijP '(e)/p = s^e = v)i(v - v)j/°(v)dv 

Thus, the course-grained kinetic equation, correctly describing isotropic and homoge- 
neous turbulence (zero spatial gradients) and satisfying the dynamic constraints pre- 
sented above, is : 
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^ + V • v/(v) + . (Vp ')/(v) + (vV^) -^ ^ /(V ) = 4 i (5.11) 

The equation (5.11) is the main result of this paper. 

5.2. Hydrodynamic approximation . 

Integrating (5.11) over v gives continuity equation (2.3). In addition, multiplying it by 
v and integrating leads to: 



(5.12) 

or 



(5.13) 

In the zeroth order, averaged with with the pdf (5.7), the equation (5.14) gives the Euler 
equation: 

dV 1 

— +V-VV = — Vp (5.14) 
at p 

where p = p(0 + JC). It follows from (5.7) and (5.11) that in the case of strongly nonequi- 
librium flow as an unperturbed state, the equation of state p'(9, p) in (5.11) is decoupled 
from the "temperature" JC entering the problem as a parameter in the zero -order pdf 
(5.7). Indeed, the coarse grained equation (5.11) gives for ^(v — V) 2 = JC : 

dpJC 



, n + ViO>V-£) + -p( Vi - Vi){vj - Vj)Sij + -V tP ( Vl - Vi)( Vj - V 3 f + - P £ 



+ l '"- V »^- V »' .V P - = 0(5,5) 
In the zeroth order (77 = 0), all gradients equal to zero, the equation (5.15) reads: 

We also define a vector 



Ei = («i-Vi)(v-V) 2 



6. Expansion. 

To find the solution to the kinetic equation (5.11) we following Chen et al (2004) write: 



and 



/(v) = /° + e/( 1 )+e 2 / (2) + - 



(6.1) 
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^- e ^ +e v^ + --- (6 - 2) 

and V = eVi. The mean of any flow property is then 

A = A^+A^+A^ + -- (6.3) 



where 



A (n) 



= J A{v).f {n) dv (6.4) 



The zero-order pdf is given by expression (5.7). The first and second -order correc- 
tions are calculated readily with i L = 3 . The result is: 

B[/(»)(v)] = [| + v V + ^ • (V P ') + (V " V ^" rf/C ^]/ (0) ^ (6-5) 
and 

- + ^ - £ ♦ v v + ^ . W) + <^2^>> - ^ 

This gives: 

(v - V) 2 /C 5/C f° 

-C- ir L -d)-V-V-2(v- V) a —} + L.(( V - V) ■ Vp - pV ■ V) (6.7) 

and 

f W f (°) (5 fl Cv - V) 2 dK. 

~ — = uc [2{v - v)a{v - v)ti{SaP - ^ Spp) + { ic - d - 2 ^ v - V ^Q^} 

f)ft (v — V"i 2 K~ f° 
-2{v -V) a £-- ( { K ] d)-X7 ■ V] + ^((v - V) ■ Vp - pV ■ V) (6.8) 

6.1. First order. Momentum equation. 
Using the pdf given by (6.8), we can write: 

ay = + 4 X) = p( k + J{6))&ij - 2 P t t IC(S 13 - 5 -fS pp ) - r T p/CV • VSy (6.9) 
and the normal stress (pressure) p = K +p'. Setting p9 = const gives; 



E = S° + S 1 =-(d + 2)^V(p/C) (6.10) 



and: 



^ + V&iVj + -V jP - (d + 2)ZdLv jP K = VjlCvrpSij + V • V) (6.11) 
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or introducing turbulence kinetic energy K = |/C: 



^ + ViViVj + -V jP -(d+ 2)^V jP K = hjjKTrpSij - <%V • V) (6.12) 
where p = p(9 + JC). 

Two remarks are in order. Firsrt, an additional dimensionless expansion parameter, i?„ = 
tt£/IC, which is not present in the ordinary Chapman-Enskog expansion, appears in the 
present noncquilibrium theory. In strong turbulence, where molecular relaxation can be 
neglected, the characteristic time is tt ~ JC/£ and i?* = 0(1) = const, is nothing but the 
renormalized Reynolds number i?» = v' rms LjvT = 0(1), corresponding to the fixed point 
of isotropic and homogeneous turbulence. Second, the derived pressure p — p(IC + 9) leads 
to rcnormalization of the sound speed by the turbulent velocity fluctuations, derived in 
Staroselsky et al (1990) and predicted earlier by Candrasekhar (1995) in the context of 
gravitational collapse of interstellar gas. 

6.2. First order. Energy equation. 
Combining the above results we, writing for the sake of simplicity of notation Sij — 
~j~Spp — Sij, have: 

^ + V • VJC = - d T T KS% -- d £ + • W + ^2 V • /CttV/C + • V (6.13) 

In what follows the last term in (6.15), responsible for the so called dilatation effects in 
compressible flows, will be neglected. Multiplying this equation by d/2 we summarize the 
above results: In the first order of the Chapman -Enskog expansion, the coarse grained 
kinetic equation (5.11) gives : 

^ + ViViVj + -V ]P - (d + rfj^VjpK = 

^VjKTT^Sij + <%V • V) (6.14) 

and 

^ + V-VK = - d r T KS% - £ + \yV . Vp > + ^t^V • KttVK (6.15) 



7. Second Order. 

According to (6.6) 

/ (2) =TTBrTjB/ (o_ rT ^!) (7-1) 

where the operator B is defined in (6.5). The symbol d tl stands for the time-derivative 
in first-order contributions to the hydrodynamic equations. For example: 
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^ = (d + 2 )^ V ^ + V,/Cr T (2^ + %V • V) (7.3) 



and 



It is clear that is a time derivative of all 0(t t ) contributions to expansion and, 
according to (7.1), = 0(r|). Thus, 

3/(0) /(Q) (v-V) 2 4 (rf + 2) ^ (rf + 2)r T g 

-fltT = ^ [( ^^~ " d ^d TTlcs v + ' Tt/cv/c P dK {Vp ] ■ (v/c) " 

2(*i ~ V;) (2(d + 2 )g|vip/C + VjJCvrpSij + V 3 - • V))] (7.5) 
Now, assuming p = const, 9 = const and V • V = 0, we rewrite (6.8) as 



' 2/C 



where 



G = [2(v V) a (v V) p (S a[i - 5j fS pp ) + ( (V -d-2)(v- V) a ^-} 



and 



f(2) R R f(0) 9 / (0) / (0) 2^2^ ,(0), ^ ^ V7V^G, ^ (0) 

Thus: 

BfW = -BvrBfM = 

(v-V) 2 <9£ o 

-^{tt[2( V - V) a (t, - + -d-2){v- V) a ^-} 2 + 

2Kt t [-2[{v - V) a (v - V)p) -(v- V) V a p - (v - V) a V pP }S a .p) + —(v- V) a (v - V)p(VT T S af3 )}} 

Tt 

-/ (0) ^(v-V)-V/C (7.6) 

where (DA) = + v • V)A Collecting all contributions up to <3((V/C) 2 ) simple calcu- 
lation with & = (v — V)i gives: 

tt ^ 



£ 2 1 <9f(°' 

(Sa,/jV 7 /C)^^^(|r - d - 1) + — &»UZ>TrSa,/3)} " "ST" + 0((V/C) 2 ) (7.7) 



so that for i ^ j: 
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2 JC 2 dV n „ BV r 



4 ?) = ^{8S W S„ - ^S a ,i - 2^S ad + "-{Vr T S t3 )} + 0(V/C) (7.8) 
and finally: 

J* _ ^pf« + %L%*. + %L%± + -IpfcrfWI + (V£> (7.9) 
The remaining 0(V/C) contributions can be easily calculated. For example : 



v(2) ttte&i? a Y\\ t t £ o vr 2(rf + 3)(rf + 5) r|g c 



We see that 

Sf } ^ 2(rf + 3)(d + 5) 



ttI^I (7.11) 



S (i)'~ d(d + 2) 

The usually neglected 0(WK) contributions to the turbulent models, may not be too 
small in the high-gradient regions of the flow. 



8. Turbulence Models Generator. 

Introducing effective viscosity 

u(]C,£, \Sij\,x,t) = vq + v t 

with 

vt = -J T ^-]C = \t t K (8.1) 
a 2 a 

gives: 

, )M , r „ _ „ _ _ . dv 2 , n dVi dVj dVj dV a dVi dV a , 

(8.2) 

We can summarize the results of this section: in the second order of expansion if p9 = 
const and V • V = 0, the turbulence model reads: 

^ + V&iVj + -V jP - (d + 2)^V jP tf = 
2 Vi ^, + v..? - 2(rf + 3 f + 5) ^^V^ (8.3) 

and 

f +V-VA' . E<".Vp- + S!^) V.^V K - V, 2 '" + f + 5 » f &J V, 

(8.4) 

In the high Reynolds number regions of the flow, i?„ = tt£ /if = const and 
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(8.5) 

The expressions (8.4) and (8.3) have a few features not present in familiar turbulence 
models. First, consider a simplest possible example of the fully developed statistically 
steady flow in an infinite channel. In this case, = 0, dtu = d x u = and the equation 
(8.3) reads: 

ldp_d_ du (rf + 3)(d + 5)r|g du OK 
p dx dy dy d K, dy dy 

The last contribution to the right side of equation (8.6) describing interaction of the 
energy flux in the wave-number space and momentum flux in the physical space, is not 
small in the proximity to the walls. This effect does not appear in familiar turbulence 
models. 

The equations (8.1) - (8.5) are quite involved and, as in the theories based on renormalized 
Wyld's expansions of the Navier-Stokes equations, accounting for the higher-order contri- 
butions to the hydrodynamic approximation is impossible. However, the main difference 
between the approach developed in this paper and the theories based on the renormalized 
perturbation expansions of the Navier-Stokes equations is that here the problem of re- 
summaion of an infinite series does arise: the equations (8.1)-(8.5) and all high-order 
models, we even cannot write down, are contained in a simple kinetic equation (5.11). 



9. Discussions and Conclusions 

The Boltzmann-equation (1.1)-(1.2) contains the Navier-Stokes equations capable of ac- 
curate description of turbulent flows. As was demonstrated in the nineties (Yakhot et al 
1992), Rubinstein (1990), ) the coarse -graining or small-scale elimination, applied to the 
NS equations, leads to the low-order turbulent models. We believe that, in principle, the 
coarse - grained Boltzmann equation (1.1) -(1.2) contains coarse-grained hydrodynamic 
approximations and all possible turbulence models (both LES and transport) of an ar- 
bitrary nonlinearity and complexity. At the present time, this coarse-graining procedure 
does not exist. One can attempt to develop the Wyld diagrammatic expansion applied 
directly on the non-linear equation (1.1)-(1.2), eliminate fast modes and, as a result, 
derive turbulence models. This well-defined program is extremely complex and we were 
unable to achieve much progress. 

Instead, the theory proposed in this paper is based on a simple qualitative relaxation 
time approximation (1.1), (1.3) combined with the 'turbulent relaxation time ' tt ~ JC/£ 
which is an exact consequence of the energy balance (4.6). Since kinetic energy is the 
large-scale property of turbulence, this relaxation time is the largest time-scale in a flow 
leading to the largest 'turbulent viscosity' vt = 0(v' rms L) = 0(JC 2 /£. This way, the fast 
modes are overdamped and the resulting equations correspond to transport turbulent 
modeling describing velocity fluctuations on the scales r w L. 

1. Transport modeling vs LES. The choice of the often - used Maxwell-Gibbs distribution 
function (2.1) is an important step in derivation of the Navier-Stokes equations from 
kinetic theory. Given these equations, one defines a spatial cut- off A and, averaging over 
the small-scale ( r < A) fluctuations, derives the coarse-grained NS equations, called 
turbulent models (Yakhot (1992), Smith et al (1992), Rubistein et al (1990), Yoshizawa 
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(1987)). By construction, the resulting equations for the "resolved scales" involve spatial 
derivatives evaluated on the mesh (cut -off), i.e. all gradients in these equations are: 



As was shown above, the high-order modeling involves high powers of velocity gradients 
(velocity differences), which, due to intermittency ( A <C L), cannot be expressed in terms 
of small deviations from the gaussian pdf. Therefore, to describe a flow, including the 
inertial range ( r s=s A <C L ) dynamics, the model must satisfy an infinite number of non- 
trivial dynamic constraints involving, among other terms, pressure-velocity correlations 
(Yakhot and Sreenivasan (2006)). With increase of the difference L — A, more and more 
constraints become relevant and, at the scales r <C L, the quality of all existing LES 
models rapidly deteriorate: with increase of the moment order, the difference between 
full (DNS) and model (LES) simulations of high-order moments grows. 
As A — ► L, the probability density of velocity difference is extremely close to Gaussian, 
the moments (u(x + A/2) - u(x - A/2) 2 ™ +1 -> 0, (u(x + A/2) - u(x - A/2) 2n -» const 
and the problem of anomalous scaling does not appear. (See Fig.). This huge simplifi- 
cation was the main and only reason for our choice of the non-equilibrium probability 
density (3.1) leading to the gaussian pdf (5.7) of the coarse-grained hydrodynamic field. 
Thus, the procedure developed here is valid for transport modeling only. To work out a 
similar approach to the LES (A C £), one has to use the expressions for the multifractal 
pdfs of velocity difference, which is not a simple task. 

2. The models. At the large scales r ss L, the zero-order (n = 0) flow can be described 
in terms of only two parameters: kinetic energy K and mean dissipation rate £, so that 
tt ~ K/£. Despite superficial similarity of relaxation times, the coarse-grained kinetic 
equation, derived in this paper, is very different from the well -known K — £ model, 
widely used in engineering: 

1. No explicitly written momentum equation is involved in our model. 

2. No separate equation for kinetic energy K is needed: all information is contained in 
the kinetic equation; 

3. The simple equation (5.11) contains K — £, Reynolds stress and all non- linear models 
etc. Thus, for example, it can describe the secondary vortices in a square duct flow, which 
no K — £ model is capable of doing. 

4. No modeling of pressure-velocity correlations is needed: it is hidden in (5.11). This 
fact can be of extreme importance for computing complex flows. 

5. The rapid distortion and memory effects are accounted for in (5.11). 

6. Sound speed renormalization. The kinetic equation (5.11) contains an equation of 
state p'(p, 9), which is a remnant of the integrated out microscopic, close-to- equilibrium, 
modes. The equilibrium speed of sound is defined as usulal: c 2 = (^)s- As follows from 
(6.12), the speed of sound obtained from hydrodynamic approximation, derived in the 
lowest order of the non-equilibrium CE expansion, proposed in this work, is rcnormalized: 
pel = 9 + K. This effect, appearing in the higher order of the renormalization group 
procedure applied to the NS equation, was derived and numerically tested in Staroselsky 
et al (1990) and, even earlier, was suggested by Chandrasekhar (1955) who was interested 
in the role of turbulence in dynamics of the interstellar gas collapse. This feature may 
prove important for simulations of turbulent flows. 

During last two hundred years, the Navier-Stokes equations enjoyed remarkable success 
in describing low- Knudsen and low- Weisenberg number fluid flows. This success was 



du _ u(x + A/2) - u(x - A/2) 
dx ~ 2A 



(9.1) 
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habitually attributed to the large scale- separation between hydrodynamic and micro- 
scopic modes: A/L < 1. There is one caveat in this argument, though. One must bear 
in mind that the NS equations are closed by imposing proper equations of state relating 
pressure to density and temperature. Invariably, the equilibrium equations of state or 
equilibrium thermodynamic relations are chosen for this purpose. The problem is that 
the CE procedure, leading to the NS equations, is an expansion around the zero-spatial- 
gradient equilibrium state of fluid and both the pressure gradient and the gradient of 
turbulent kinetic energy are large -scale properties defined on the same hydrodynamic 
scales. Thus, in this state of fluid, the only relevant dimensionless parameter is Mach 
number Ma = v rms /c s . Therefore, to describe the stresses, in addition to thermody- 
namic considerations based on an estimate p = 0{pc 2 ), one has to add the turbulent 
contribution 8p w pv 2 ms , which appears naturally in the procedure developed above. 
Thus, strictly speaking, the pressure term in the NS equation, which is the functional of 
the solution, must be determined self-consistently. 

Given the relaxation time tt ~ K/E and effective viscosity vt — ttK, the kinetic 
equation (5.11) generates various fields such as V and K. However, to evaluate the 
relaxation time, the magnitude of the local dissipation rate is needed. This can be done 
selfconsistently, if, in accord with Kolmogorv theory, we assume 



e = M y^r (9-2) 

Keeping in mind the Lattice Boltzmann applications, we can define the velocity deriva- 
tives on a lattice of lattice constant A with the result written for simplicity in a one- 
dimensional case: 



g(Q M K Vm " V< ^ V(i) ) a (9-3) 

where i denotes position (coordinate) of a lattice site and the averaging is carried out 
selfconsistently on a non-equilibrium pdf, which is a solution to the kinetic equation 
(5.11). Another possibility is to use the dissipation rate directly from the energy balance 
(8.5). 

The proportionality constants in (8. 2) are reasonably close to those quoted in Speziale 
(1987), Rubinstein ct al (1990) and Yoshizawa (1987). It should be pointed out that the 
value of calculations based the low-order trancations is not clear. Indeed, recasting the 
results of Sections 7 and 8 in terms of dimensionless rate of strain r\ w K\Sij\/£ , we find 
that cr^ 2 ) w v 2 S 2 /K w Kr) 2 and the next order gives: cx Krf . Taking into account 
that in the flows of engineering importance the parameter rj is not small, often reaching 
values 77 « 10 — 20, the perturbation expansion does not converge and no trancation is, 
in general, possible. We would like to reiterate that the entire series is contained in a 
relatively simple kinetic equation (5.11). 

The low-order hydrodynamic approximations and turbulence models, contained in (5.11), 
are similar to the well-known ones, which have been widely tested in both scientific 
and engineering environments (see Chen (2004), for example). The future simulations of 
strongly nonlinear, rapidly distorted or oscillating flows, will provide a decisive test of 
the limits of validity and accuracy of calculations based on the kinetic equation derived 
in this paper. 

I am grateful to H. Chen, A. Polyakov, I. Staroselsky, X. Shan, S. Succi, E. Kamen- 
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skaya, T. Gatski, K.R. Sreenivasan and U. Frisch for most interesting and stimulating 
discussions. 
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